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Abstract 

We derive a generalization of the Destri - De Vega equation governing the scaling 
functions of some excited states in the Sine-Gordon theory. In particular configurations 
with an even number of holes and no strings are analyzed and their UV limits found 



to match some of the conformal dimensions of the corresponding compactified massless 
free boson. Quantum group reduction allows to interpret some of our results as scal- 
ing functions of excited states of Restricted Sine-Gordon theory, i.e. minimal models 
perturbed by (friz in their massive regime. In particular we are able to reconstruct the 
scaling functions of the off-critical deformations of all the scalar primary states on the 
diagonal of the Kac-table. 



1 Introduction 



The exact computation of scaling functions in two-dimensional integrable quantum field the- 
ories has been paid a lot of attention in the last years. The Thermodynamic Bethe Ansatz 
(TBA) methods JT|, H, for example, are able to provide, for a theory whose factorizable S- 
matrix is known, a very good framework to compute a set of coupled non-linear integral 
equations driving the evolution along the Renormalization Group flow of the Casimir energy 
of the vacuum on a cylinder. Towards the UV limit this quantity is directly related to the 
central charge of the corresponding Conformal Field Theory. 

However, the TBA method suffers of some disadvantages. First of all the deduction of 
the integral equations is done starting from the Factorized Scattering Theory (FST) and not 
from the Quantum Field Theory (QFT) itself. Now, the FST is often only conjectured and 
the TBA itself is therefore not fully deduced from the QFT. Another point is that there are 
models where the number of equations for the TBA becomes infinite or increases to a point 
that numerical integration becomes untractable. 

An undoubtable progress of these exact methods came with the introduction, at least for 
the Sine-Gordon model, of the so called Destri - De Vega (DdV) equation f||, |J, whose major 
advantages are: 

• to summarize in one complex integral equation a large number, sometimes an infinity, 
of real coupled integral equations of TBA type; 

• to be deduced from a lattice regularization of the model itself and not from the scat- 
tering theory corresponding to it. 

Unfortunately the DdV equation is developed and well explored only for the vacuum of 
the Sine-Gordon model and it would be very important to write down analogous equations 
for other integrable field theories too, as well as for excited states. In this letter we make 
progress towards the study of DdV equation for excited states in the Sine-Gordon model. The 
quantum group reduction allows then to extend these results to the Restricted Sine-Gordon 
models, i.e. the integrable theories obtained by perturbing the conformal minimal models by 
their relevant operator $13. This is the first case of extension of the DdV equation out of the 
pure Sine-Gordon realm. 

In sect. 2, after recalling some general properties of the SG model on a cylinder, we 
write down, following |5[], the DdV equation for vacua with twisted boundary conditions. 
These vacua are states that can be obtained from the true vacuum by application of spin 
fields. In sect. 3 we observe that after quantum group reduction, some of these states can 
be reinterpreted in the framework of RSG theory as the off-critical deformations of scalar 
primary fields along the main diagonal of the Kac-table. A very precise numerical check 
against known cases confirms this procedure. In sect. 4 we afford the problem of writing the 
DdV equation for excited states over the twisted vacua, by exciting Bethe strings and creating 
holes over the Dirac sea of real solutions. This corresponds to approaching the problem in a 



true Bethe ansatz spirit. In this short letter we shall make the choice to illustrate for technical 
simplicity only the case of pure hole excitations, leaving the more general problem of string 
excitations to a future, more extensive publication ||. In sect. 5 we compute the energy as a 
functional of the solution of the DdV equation. In sect. 6 we investigate the ultraviolet (UV) 
limit of the energy and compare with the operator content of the compactified free boson 
CFT that represents the UV limit of SG theory. Some preliminary numerical investigations 
are reported in sect. 7 and in sect. 8 we draw our conclusions and perspectives for future work. 



2 The Sine-Gordon Hilbert Space on a cylinder 

The Sine-Gordon (SG) model, with action 

S = j d 2 x{\: (8(f)) 2 : +A : cos/50 :} (1) 



is invariant under the symmetry 0^0 + 27m/ (3. States transforming as — > e ia \^) 
under this symmetry constitute a so-called a-sector H a of the Hilbert Space of states 7i = 
®--K<a<-K^a- On a cylinder of radius R the degeneracy among the a- vacua is removed and the 
various a-sectors are selected by imposing suitable twisted boundary conditions on the space 
direction of the cylinder. For p = g7r _^ 2 rational, the values of a are quantized; if p is integer, 
i.e. the case that we shall mainly consider in the following, then a = k — 0, 1, ...,p. The 
a-sectors have of course a state of lowest energy, the so called a-vacuum. For a = this is 
the true vacuum of the theory, the other a-vacua can be excited from the true vacuum by 
particular operators which are called spin operators. In the CFT representing the ultraviolet 
(UV) limit of SG at a fixed p (SG P ) (a c = 1 massless free boson compactified on a circle 
of radius these operators are represented by vertex operators of conformal dimensions 
A « = 4 P (p+i) • For more details see, e.g., ref. 0. 

A DdV equation can easily be written for all the a-vacua along the lines of the standard 
deduction of the original work of Destri and De Vega |3], |j] . Indeed, one simply has to consider 
twisted boundary conditions, for which a DdV equation is written in |4j for the XXZ chain. 
It is a straightforward exercise to rewrite it in the form more appropriate for the SG P model. 
The resulting equation is P 

/+0O 
d&G{6 - ff)$sm log(l + e fa{e+l0) ) (2) 
-oo 

where r = mR, m being the soliton mass. The kernel G{6) = —^-^\ogS(9) is determined 
by the soliton-soliton scattering amplitude 

i r+°° ,,sm(9k) sinh^(p-l) 



S(9)=exp- dk — )—!- i- (3) 

K ' P 2 J-oo k sinh cosh ^ V ; 



Once the pseudoenergy f a {6) is determined by solving the DdV equation, it can be used to 
compute the ground state Casimir energy E a {R) = — , where 

c a (r) = — / ^sinh^mlog(l + e /a(e+i0) ) (4) 

7T Z J-oo 



The adimensional function c a (r) is known as scaling function for the a-sector ground state. 
A standard dilogarithm calculation shows that the UV limit r — > of this function is given 
by 

Ca(0) = c uv ~ 12(A+ - A_) = 1 - (-)* (5) 

P+l \7T/ 

thus reproducing the correct central charge c = 1 of the free boson and the correct conformal 
dimensions A + and A_ of the spin operators mentioned before. 

One can thus follow the renormalization group destiny of the a-vacua in the different 
sectors on the cylinder in an exact way, overcoming the need of approximante methods as 
the Truncated Conformal Space Approach (TCSA) [|J. The obvious result that one can see 
from a plot of the functions c a (r) for all a's is that they all accumulate towards infrared to 
the vacuum of the theory. 



3 RSGp spin states 

An interesting application of this result is the possibility to reinterpret it in the light of the 



quantum group truncation occourring at rational p, which, as shown by |10[ , selects inside 
the Hilbert space of the full SG P model, a certain subspace with a consistent subalgebra 
of operators acting on it in such a way that it can be thought as a sort of smaller theory 
embedded in the larger SG P one. The SG vacuum never appears as a state in this restricted 
Sine-Gordon (RSG P ) model. It is known that the role of the vacuum is instead played in this 
context by the a-vacuum with k = 1 (we restrict here for simplicity to the integer p case, 
it is not difficult to generalize to any rational p). Indeed, lim r ^ c 7r/p( r ) = 1 — p(p+i) ' as ^ 
must be for the correct scaling function of the vacuum of the corresponding minimal model 
perturbed by 01^. To support this result we made numerical comparison of the scaling 
function obtained by this DdV equation with that obtained from the traditional TBA of 



ref. ||11|| . The two functions agree with 14 significant digits, for integer values of p ranging 
from 3 to 12. 

The quantum group truncation eliminates only two a-vacua: the k = 0, 7r ones. All the 
other a-vacua have an interpretation as RSG P states, namely they correspond to the (off- 
critical evolution of) the primary states on the diagonal of the Kac-table of the p-th minimal 
model, the so called spin states. Indeed the computation of the UV limit of the scaling 
function gives in general c K7r / p (0) = 1 — p7~pij — 24 ^^^ , thus providing the conformal 
dimensions of the states in the diagonal of the Kac-table. The scaling functions of a few 



of these states were already explored by other authors with other means [|12|, [13|, TJJ. We 
succeded in various cases to compare our numerical results with theirs, getting the same kind 
of good agreement up to 14 significant digits as for the vacuum. 

The quantum group reduction for other states, especially for secondaries, is more involved 
and a full treatment goes out of the scope of the present letter. Let us only comment here 
that unfortunately none of the hole excitations we are going to describe in the following is in 



the class (type II representations) that survives quantum group reduction. 



4 DdV for pure hole excited states 

In this section we present our derivation of the DdV equation describing pure hole excitations 
obtained from the vacuum by creating some holes in the distribution of real roots. 

Let us consider the Bethe equations for the SG P model discretized as a 6-vertex model on 
a light-cone lattice with 2N sites. Also, take the corresponding inhomogeneous spin chain 
have tu-twisted boundary conditions, related to the previous a by a = —u^—. These Bethe 
equations have the form 

" sinMA, - X m + ij) = _ 2 JZ mih(Aj + ie B + $) to 
ii x sinh(A i - A m - ij) ii sinh(A J + iQ n - f ) 1 ' 

Here 7 = The fundamental state (a- vacuum) is the one with M = N. The light-cone 
approach corresponds to the choice of inhomogenuities iQ n = (— l) n+1 G. Define the function 

sinh(ix + A) 

0(A, x) = 1 log — — (7) 

smh{tx — A) 

and the M real solutions counting function 

M 

Z M (X) = N[4>(\ + 9, |) + 0(A - e, i)] - x: 0(A - A fc , 7) - w (8) 

k=l 

Then, the logarithm of the Bethe equations can be written asQ 

Z M (\ j )=2TrI j , IjeZ + l , j = l,...,M (9) 

We consider here a configuration with M < N real solutions and if holes (i.e. H values of 
Ij are missing). Consistency of the Bethe equations requires that the number H must be 
even [[T5|, |TB| . The range of possible values of Ij actually increases due to this "perturbation" 
of the system. Now there are in total T = M + H = N+ y possibilities, while the roots are 
M = N — y. In order to implement the usual DdV trick, we have to rewrite the counting 
function Zj^{X) by adding and subtracting the contribution of holes 

T H 

Z M {\) = NQ{X) - uj - ]T 0(A - A fc , 7) + E <K A - 7) (10) 

k=l h=l 

where we have indicated for short Q(X) = 0(A + G, |) + 0(A — 9, |). The standard derivation 
of DdV eqautions makes use of the trick 



In this paper we are interested only in states with M even 



to write (the integration path T contours the whole real axis counterclockwise) 

f>(A-A fc , 7 ) = /^(A-^Aiog^ + e^M) (12) 
= i{X*L){X)+i{X*Z M ){X) (13) 



k=l 



where L(A) = log(l + e iZM( - x+iri ^) -log(l + e~ iZM( - x ~ iriS> ) for i] — > 0, * denotes standard convolu- 
tion (f*g)(x) = f(x — y)g(y)dy and we have introduced the notation X(X) = ^0'(A, 7). 

The rest of the derivation repeats the same steps as the standard one for the vacuum 
DdV H U, with the only modification that NQ(X) must be now substituted with A = 
NQ — to + J2h=i<j>(^ ~ A^, 7). This goes to modify the imi^sinh^ term because instead of 
computing N(K * Q)(X) we have to compute now (K * A) (A), where K(X) is the inverse of 
(<5 + X)(A), that is 

[K * (6 + X)](X) = 6(X) (14) 

We have to compute terms like K * <p(X — //, 7). Just doing a Fourier transform we get in the 
case of (jl real: 

K * <f)(X - (1, 7) = -i log S(z(X - p)) (15) 

All these calcualtons are easily performed in Fourier transformed space, by remembering that 
the Fourier Transform of <j)(X, 7) is given by (j>(k, 7) = sign( 7 )|^^^^^. 

The continuum form of DdV equation is recovered if we let 9 = A7r/ 7 , a = and 
f(6) = ]im N _ f00 iZ M (X) 

H 

f(6) = imR sinh 6 + ia + {6 -8 h ) 

h=l 

/+00 
d9'G(9 - 6')Sm log(l + e /(e (16) 
-00 

where 9 h = ^X h 

The condition we use to determine the parameters 9h is the fact that they are missing roots, 
i.e. although they are not taken as roots of the Bethe equations, nevertheless they satisfy 
Zm(Xh) = 2irlh for some chosen Ih (the set of Ih for h = 1,...,H then labels the excited 
state we are examining). In the continuum limit this condition reads as f(0h) — 2mlh. To 
numerically integrate the DdV equation we start from the vacuum evaluation of f(9) and 
search for values in the 9 complex plane where f{6) = 2irlh- We then insert this rough 
determination of 9^s in the DdV equation, and iterate to get a better determination of /. 
With this better determination we recompute the values of 9h and then back to the DdV for 
a newer determination of / and so on iteratively. Although the convergence of this procedure 
slows down if compared to that of the vacuum case, it is still possible to determine the 
functions without much effort to some 7 or 8 significant digits. 



5 The energy 



The energy of the state as a functional of / can be obtained as the logarithm of the eigenvalue 
of the diagonal to diagonal transfer matrix corresponding to the state at exam 

£m = ^X>(^)-2tt] (17) 

where $(A) = 0(A + 6, ^) — 0(A — 6, ^). As before, we add and subtract a contribution from 
the holes 

N T H 
E M = -5 E - E $ (^) - 27rAf] (18) 

The first sum in this expression can be manpulated in a manner similar to what done in the 
deduction of DdV equation 



T 



X>(A 3 ) = ^* M _log(l + e < 

/+oo /7<r» /*+oo /7<t> 

— $'(x)Z m (x) - i / — $'(ar + iry) log(l + e^^-H")) 
-oo 27T J-oo 27T 



-oo 

-°° do: 



+ i / — $'(x - ir]) log(l + e^*""") (19) 

J-oo Z7T 

Now we substitute the Zm{^) appearing in the first integral of this expression using the DdV 
equation (|16|). Also, rewrite the last term 2ttM as 2tt(N — -j). Introducing the notations 

V 2 ( n _ f+oo ; ^(A + 26,f) > 



£, = ¥(" 2l + 7-oo rfA ^ShjA ' (20) 

and 

■0(A) = ^[sech(^(A - 6)) - sech(^(A + 6))] (21) 



we can now write the energy as follows 

/+oo 
dAV>(A + i77)log(l + e 
-oo 



+ % / d\^{\ - irj) log(l + e -«M*-»»)) 

N JL r+°° dx N N JL 

+ ijj E log " A*)) + ^ - - E <&(A fc ) (22) 



The first term E c reproduces the standard bulk energy |3J. The second term, in the N — ► 
oo limit acquires the same form as the vacuum energy in ref. [Q. Notice, however, that 
this similarity is only formal, because now this term is evaluated with the "excited" f(9). 
The physical interpretation is that of an energy of the vacuum distorted by the presence of 
particle excitations. The last three terms are new and, when evaluated in the N — > oo, = 
-log — — > oo limit, concur to give, after some algebra, a term of the kind mELcosh^, 



with the clear interpretation as a particle energy. Concluding, the continuum limit energy 
E(R) of the excited state, epurated from the bulk energy and the infinities which are artifacts 
of the lattice, is 

H fji r+oo 

E(R) = lim (E M - E c ) = m V cosh 9 h / d9 sinh 9%m log(l + e f{e+l0) ) (23) 

Note that in the infrared limit r — > oo, where the convolution integral can be dropped, an 
if -hole state can be interpreted as an if-particle state. Having in our approach H even only, 
we see states with an even number of particles only. To deal with the more general problem of 
any number of particles, one should also consider DdV equations which are continuum limits 
over lattices with odd number of sites, or equivalently choose twisted boundary conditions in 
an appropriate way. This is, however, out of the scope of the present letter. 

6 Kink DdV equations and analytic ultraviolet limit 

The scaling function c(r) = —QRE(R)/ir is particularly interesting because, in the ultraviolet 
(UV) limit r — > 0, it provides the value of the conformal dimensions A + + A_ and therefore 
allows the identification of the states with those of the CFT Hilbert space. The calculation 
of this limit can be done in an analytic way by resorting to the so called kink DdV equations. 
Let x = log(r/2). The UV limit is then given by x — > — oo. The shape of the function —if (9) 
in the region x <C 0, is almost flat (and its value very close to a) inside the region \9\ <C —x 
while it has a sharp e rd behaviour for 9 3> — x or — e rd for < i. The hole positions 9h 
are given by the condition —if(9h) = 2nl h . F° r x — ¥ — 00 therefore, 6 h — > sign(J /l )oo. Let 
us divide the set of holes into those with positive rapidity (//, > 0) and those with negative 
rapidity (//, < 0). Call these two sets X + and X_ respectively. We assume for simplicity that 
these two sets contain half of the holes each, \T±\ = H/2. (Remember that if the number of 
sites of the lattice is even, the number of holes also must be even). The most general case will 
be treated in 0. It is obvious that towards the UV limit, the holes in X + (resp. X_) behave 
as right (resp. left) kinks. Their rapidities can be shifted as follows: ph — 9h =t x if h E T±. 
The DdV equation splits in the UV limit into two right (+) and left (— ) kink DdV equations, 
correspondingly to the general shift 9 ± x = p. Dropping terms exponentially small when 
x — > — oo and introducing the two kink pseudoenergies f±(p) = lmx^.oo /(±(p — x)), these 
two equations turn out to be scale invariant, i.e. indipendent of x 

T if±(p) = e ±p ± <*+ f a' T i £ logS(p - p h ) 

hei± 

/ + 0O 
dp'G(p - pT)Qm log(l + e f±{p+i0) ) (24) 
-oo 



where we have used the fact that log £(±00) = ±ia' with a' = rt 2 ^-. Correspondingly, the 
scaling function assumes the form c(0) = c + + c_ , with 



6 

c± = 

7T 



1 r+00 

e ±ph =f - / dpe ±p Qm log(l + e 
hex± 71 J -°° 



/±(p+iO)N 



(25) 



Introducing ip±(p) = e ±p ± a + |a' =F « 2/iez± l°g ^(P — P/i) we can write eq. (|25"D as 
c± = --Y j e ±ph + -^ dpv'±{p)%m\og{l + e s±{p+l0) ) 

IT hex± vH J-00 

P. r+00 A 

± "2 E / ^— log^p-p^mlog^ + e^^ )) (26) 

The second term on the r.h.s. of this equation can be treated by resorting to the lemma 
cited in ref. [§] and observing that /±(=Foo + iO) = ±za^-. The first and third terms can 
be instead manipulated substituting e ±ph by use of the DdV equation and remebering that 
Ph — Ph' — > for x — > —00 if both h, h' are in X + or X_ and log 5(0) = 0. One arrives at the 
final formula 

c(0) = 1 - 12(A+ + A_) (27) 

with 

(f(p+l)T«) 2 , „ „ +v , ff 2 , 2R , 



giving the left and right dimensions of the conformal state which is the UV limit of the H hole 

Hk 
2p 



state we are examining. The (Lorentz) spin of the state is given by A + — A_ = — ^r+J2h=i ^h- 



The numbers K± are non-negative integers, and account for secondaries. 

The states we are studying are uniquely defined by giving the sequence of I2, Ih)- 
Repetitions of values of Ih must be avoided in a sequence. As an example, we discuss here the 
case p = 3. The k = sector corresponds, at the UV limit, to conformal states selected by 
periodic boundary conditions. They are encoded in the modular invariant partition function 
of the model. An inspection of this latter shows that it contains a certain number of states 
of spin 0, among which there are the A + = A_ = |,|,3 states. These can be obtained in 
our framework by exciting 2,4,6 holes respectively, with the lowest possible values of \Ih\ 



1 1\ 1 (_3 _1 1 3\ 4 (_5 _3 _1 1 3 5n 

2' 2/ 3 ' V 2' 21 21 11 3 ' V 2' 2' 2>2>2>2' 



(29) 



Using larger values of one can get secondaries of these states, even with spin different 
from zero. Often, there are many degenerate secondaries with the same values of A + + A_. 
We can also obtain the same values of A + + A_ by different sequences of J^'s. In the example 
above 2- hole states and 4- hole states often degenerate in the conformal limit, although their 
infrared destiny is very different, going to 2-particle and 4-particle states respectively. To link 
the description of this degenerations in terms of hole states with the Virasoro (or extended 
conformal algebra) description is a very interesting problem to investigate. 



Of course the modular invariant partition function mentioned before, or its analogs for 
other values of p, contain other states of spin that we cannot reach with our construction as 
far as we do not deal with string excitations (see |§). The hope is, once both hole and string 
states are taken into account, to show that the whole modular invariant partition function 
can be reconstructed. 

If one considers now hole excitations over the other a-vacua with k ^ 0, one should expect 
to reconstruct, in the UV limit, the partition functions with a-twisted boundary conditions. 
It is a straightforward exercise of CFT to write down these partition functions and compare 
the operator content of these sectors with our results. The states we find by choosing 2,4,6... 
holes with the "minimal" choice of I^s have now fractional spins. It is a matter of fact that 
such exotic states are indeed present in the twisted partition functions. 

7 Numerical work 

The results of the previous section have been checked by numerical integration of the DdV 
equations. Apart the already mentioned numerical study of the a-vacua, that confirms with 
high precision the expected identification with RSG P excited states on the diagonal of the 
Kac-table - a new result by itself - we have investigated numerically the hole excitations with 
2,4,6 holes for values of p = 3, 4, 5. The precision, in this investigation has been taken a bit 
lower, say 6 or 7 significant digits. It is a simple choice of speed, no technical problem prevents 
to get higher precision data. A detailed summary of these numerical results, together with 
their particle interpretation, level crossing phenomena, etc... will be given in ||. Here we 
only present, as an example of the method, the energy evolution against R of few of the two 
and four holes spin zero excited states over the a = vacuum in p = 3. The self explanatory 
Fig.l collects these data. One can clearly see that the if- hole states accumulate to if -particle 
states in the R — > oo limit, as commented before. 

It is interesting to note the undoutable existence, in this graph, of a level crossing between 
the (— |, |) and the (— |, — |, ~, |) states at R = 5.613.... Such level crossings must be present 
in an integrable model. One of the advantages of knowing scaling functions exactly, compared 
to other approximate methods of investigation like e.g. TCSA, is that the existence of a level 
crossing can be established with no doubt. There is no danger here to confuse a level crossing 
with two lines that approach very close and then repulse each other without crossing: each 
line, in the DdV approach, is computed indipendently from the others. 

Of course there are many states missing in Fig.l. We have not investigated all the states 
that include strings in their Bethe root configurations. These should complete the picture 
and most probably provide for more level crossings. We have also done some investigation in 
the other a-sectors, that we do not reproduce here for pure reasons of space. In these sectors 
the number of spin states is much lower. It is interesting to notice the existence of fractional 
spin states. The a = ir sector is particularly interesting, because it realizes the antiperiodic 
boundary conditions where it should be possible to see the spin | states corresponding to the 




Thirring fermions. 

8 Conclusions 

We have considered here hole excitations of the SG model over the a-vacua, finding a DdV 
equation for them and calculating their exact scaling functions, numerically for generic R 
and analitically for R — * 0. In this latter limit, we have found agreement with the operator 
content dictated by the CFT partition functions in the variuos twisted sectors. Of course, we 
have reproduced only a part of the spectrum, because we have not considered more compli- 
cated states including strings. We intend to deal with the general problem in a forthcoming 
publication JJ. In this short letter we made the choice to deal with this technically simpler 
problem to probe the method and show its efficiency The only states we have provided 
for RSGp models through q-reduction are for the moment those corresponding to the alpha 
vacua, i.e. the off-critical deformations of the spin primary states which lie on the diagonal 
of the Kac-table. It is of course obvious that once the full set of SG states will be treated, 
access to the scaling functions for all the RSG P states will also be available. The extension of 
the method to non-integer rational p could then allow the treatment of non-unitary minimal 
models perturbed by 0i3 too. Among these, the Lee- Yang model plays the particular role 
of being the simplest model, on which both TCSA data and recent TBA exploration of few 
excited states are available |T7], ITBf. It would be extremely interesting to make a comparison 



between these results and our approach. In particular the method of ref. [H| suggests to 
explore monodromy issues of the scaling functions also in our approach. The link between 
monodromies and Bethe ansatz interpretation of the results would certainly provide more 
insight in the problem of excited states in general. 

Acknowledgements - We are endebted to F.Capocasa, CDestri and P.E.Dorey for 
useful discussions. G.Vacca is also thanked for interest in this work. E.Q. thanks the warm 
hospitality at Physikalisches Institut der Universitat Bonn during the early stages of this 
project. This work was supported in part by NATO Grant CRG 950751. 

References 

[1] Al.B.Zamolodchikov, Nucl. Phys. B342 (1990) 695 

[2] C.N.Yang and C.P.Yang, J.Math.Phys. 10 (1969) 1115 

[3] C. Destri and H. De Vega, Phys. Rev. Lett. 69 (1992) 2313 

[4] C. Destri and H. De Vega, Nucl. Phys. B438 (1995) 413 

[5] Al.B.Zamolodchikov, Phys. Lett. B335 (1994) 436 

[6] D.Fioravanti, A.Mariottini, E.Quattrini and F.Ravanini, in preparation 

[7] Al.B.Zamolodchikov, Nucl. Phys. B432 (1994) 427 

[8] V.P.Yurov and Al.B.Zamolodchikov, Int. J. Mod. Phys. A5 (1990) 3221 

[9] F.Smirnov, Nucl. Phys. B337 (1990) 156 
[10] D.Bernard and A.LeClair, Nucl. Phys. B340 (1990) 721 
[11] Al.B.Zamolodchikov, Nucl. Phys. B358 (1991) 497 
[12] PFendley, Nucl. Phys. B374 (1992) 533 
[13] T.Klassen and E.Melzer, Nucl. Phys. B370 (1992) 511 
[14] M.J.Martins, Phys. Rev. Lett. 67 (1991) 419 
[15] H.De Vega, Int. J. Mod. Phys. A4 (1989) 2371 
[16] M.Karowski, Nucl. Phys. B300 (1988) 473 

[17] V.V.Bazhanov, S.L.Lukyanov and A.B.Zamolodchikov, preprint CLNS 96/1416, LPM- 
96/24, |hep-th/9607099| 



P.E.Dorey and R.Tateo, Durham preprint DTP-96/29, [hep-th/9607167 



